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Abstract 

We report calculations of the ground state energies and geometries for clusters of different sizes 
(up to 80 particles), where individual particles interact simultaneously via a short-ranged attractive 

o. 

O , - modeled with a generalization of the Lennard- Jones potential — and a long-ranged repulsive 

Yukawa potential. We show that, for specific choices of the parameters of the repulsive potential, 

£T) ■ the ground state energy per particle has a minimum at a finite cluster size. For these values of 

the parameters in the thermodynamic limit, at low temperatures and small packing fractions 

where clustering is favored and cluster-cluster interactions can be neglected — thermodynamically 
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stable cluster phases can be formed. The analysis of the ground state geometries shows that the 
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spherical shape is marginally stable. In the majority of the studied cases, we find that, above a 
certain size, ground state clusters preferentially grow almost in one dimension. 



I. INTRODUCTION 

Understanding the formation of self-assembled structures at nano and mesoscopic level 
is one of the central issues in condensed matter studies |lj, and particularly in biological 
and soft matter fields. More recently, significant efforts have been made in the direction of 
finding connections between features of the inter-particle potential and the resulting stable 
supra-particles patterns. In this respect, an interesting class of systems in which particles 
self assemble into large aggregates is generated by systems of particles interacting simul- 
taneously via an attractive potential (which can be of van der Walls or depletion origin) 
and a repulsive interaction (usually of screened electrostatic origin). In density and tem- 
perature conditions where the attractive part of the potential would generate macroscopic 
phase separations, the presence of a repulsive part contrasts the phase separation process, 
leaving the system, in some cases, in a microphase separated condition. The competition 
between attraction and repulsion generates the formation of equilibrium phases where stable 
clusters of particles are formed (sometime called cluster or micellar phases) or even of more 
complex structures (such as lamellar or columnar phases). Experimental evidence for such 
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yjpresented for colloidal systems 
. Cluster phases have also 



thermodynamically stable cluster phases have been_recent 

and for solutions of proteins at low ionic strength 

been observed in aqueous solutions of silver iodide |7| . Theoretically, the formation of phases 

with different morphology has been addressed within lattice models and mean field calcu 

lations for „e d e.ect^atic interactions QQU 0QQ and, mfiE e recent., in 
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colloidal science, explicitly accounting for ion condensation around particles L 1 J • Numerical 
evidence of cluster phases has also been recently reported |l2|, |l5|, [la [l7| . Interestingly, at 
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low densities, particles aggregates can act as monomers of a supramolecular 
phenomena typical of the liquid state as crystallization and glass transition 

In this work we study an off-lattice three dimensional representative model of particles 
interacting via a short range attraction (which we model with the generalization to very large 
a of the Lennard- Jones (LJ) 2a — a potential, as previously proposed by Vliegenthart et 
al. |l8[) complemented by a long range repulsion, which we model by a screened electrostatic 
Yukawa potential. Screening is included in the model, since we are mostly interested in 
describing charged colloidal and protein solutions, where addition of salt in the dispersant 
medium provides an efficient mechanism of reduction of the repulsive interaction range. 



Modeling the short range attraction with the Lennard- Jones 2a — a potential is of course 
arbitrary, but at the same time representative (by tuning the exponent a) of all short- 
range potentials, arising either from depletion interactions or from van der Wall forces. For 
different choices of the repulsive potential parameters, we calculate the ground state energy 
and structure of clusters of different size, to evaluate the conditions under which stable 
cluster phases are expected at low temperatures, where entropic effects can be neglected. 
We show that, on varying the range of interaction and the intensity of the long range repulsive 
potential, it is possible to obtain stable cluster structures of different geometry. To locate 
the ground state structures and their energies, we have used a modified version of the basin- 



hopping algorithm 19] which has been proposed and extensively applied in models of simple 
liquids and molecular systems. To extend the range of studied cluster sizes beyond what it 
is possible with present numerical resources, we complement the ground state calculations 
by analytic calculations of the cluster self energy, under the assumption of spheric cluster 
shape. 

The paper is organized as follows. Sec. |TT1 contains a study of the clusters ground states 
for the 2a — a attractive potential for two values of a, typical of short-range attractive 
interactions. A comparison of ground state energies and cluster shapes between the short 
range attractive potential and both the standard a = 6 case (from Ref. [19() and the Hard 
Sphere (HS) case (from Ref. |20j ) is also reported. In Sec. IHII we present cluster energies and 
geometries resulting from the addition of a repulsive Yukawa potential, for different values 
of the potential parameters. Sec. HYI is devoted to the analytic calculation of the cluster 
self-energy under the assumption of spherical shape. Finally, Sec. El contains a discussion 
of the results and conclusions. 

II. CLUSTER GROUND STATES FOR ATTRACTIVE SHORT RANGE 2a - a 
POTENTIAL 

Short-range attractive interactions arise in colloidal systems either from van der Waals 
interactions (integrated over the volume of the interacting particles) or via a depletion 
mechanism, i.e., when particles of intermediate dimension between the suspended colloids 
and the solvent molecules are added to the solution |21| . Here we model the short-range 
attractive potential in a numerically convenient way with the generalization of the Lennard- 



Jones potential, namely the 2a — a potential, previously proposed by Vliegenthart et al. 
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where the value of a, which in standard LJ case is fixed to six, is varied to control the range 
of interactions. In what follows we have chosen e and a as the units of energy and length 
respectively. We focus our attention to three values of the parameter a that are relevant for 
our discussion. Firstly, a = 100, which corresponds to an attractive range of the order of 
a few percent of a, typical of vary narrow-ranged depletion interactions. Next, we consider 
the case a = 18, which corresponds to the case for which the value of the potential energy 
at the distance of the second neig hbor shell becomes negligible. Finally, for comparison, we 
also report results from Ref. |l9j], for the Lennard- Jones case a = 6. The three potentials 
(a = 100, 18, 6) are illustrated in the inset of Fig. [TJ 

To calculate the ground state energy of clusters of different size N and their geometry, 
we have used a modified version of the basin- hopping algorithm [19l| . introduced by Wales 
and Doye. Such an algorithm consists of a constant-temperature Monte-Carlo simulation 
where the actual potential energy surface is transformed into a stair-like surface by replacing 
the potential energy function with the values of the potential energy of the closest local 
minimum, named inherent structure [22]. Acceptance criterion is thus based upon the 
change in the inherent structure energy. In this way, particles can be moved to regions with 
high potential energy with large acceptance rate, facilitating the overcome of the potential 
energy barriers. To further favor barrier crossing, every 100 Monte-Carlo steps, the less 
bounded atom is removed and re-inserted in the position with lowest insertion energy. The 
largest cluster size studied is of the order of 80 particles. 

Fig. H shows the ground state energy per particle E/N as a function of iV -1 / 3 , for all 
three values of a. Data are reported as a function of iV~ 1//3 since, for large clusters, one 
expects to observe the functional form E/N = Cq + ci/iV 1 / 3 . Indeed, the cluster energy can 
be expressed as the sum of a bulk term, CqN, and a surface term, C\N 2 ^. We find that, 
for a = 100, when clusters are larger than « 13 — 15 particles, the ground state energy per 
particle is well represented by the fit E/N = — 5.937 + 7.435/iV 1 ' 3 , providing an estimate for 
the bulk energy compatible with a close-packed fee or hep crystalline ordering value (—6). 
It is interesting to note that the a = 18 data are practically indistinguishable from the 
a = 100 case, suggesting that cluster ground state energies and geometries are the same for 



all potentials with a > 18. 

Fig. ^ also reports the ground state energy per particle for the standard Lennard- Jones 
case a = 6, from Ref. |l9(. As compared to the shorter range case, the ground state energy of 
the L J is lower, due to the attractive contributions arising from second and further neighbors. 
The asymptotic behavior (with iV^ 1 / 3 ) is approached at much larger sizes (order of ~ 30 — 40 
particles) and suggests an extrapolation to the bulk value, consistent with the estimate for 
the LJ fee or hep values of —8.61 [23]. 

The sign of the surface term, i.e., the sign of C\, provides information on the thermody- 
namic behavior of a macroscopic system at low temperatures. Indeed, when T is low, the 
ground state energy becomes the relevant thermodynamic potential. The condition c\ > 
ensures that the system ground state in the thermodynamic limit is composed by a sin- 
gle macroscopic bulk cluster. In the cluster-based thermodynamic description pioneered by 



Hill 24], C\ > indicates that at low temperatures the system will phase separate in a dense 
bulk liquid in equilibrium with a gas phase. 

Next, we discuss the geometry of the ground state clusters. Fig. El contrasts the cluster 
geometries for the short-range case with the LJ case. For very small N, the structure of the 
clusters is independent on the range of the potential. For iV < 10 the structure of the clusters 
does not change with a. For larger iV values, the geometry changes significantly, since the 
constraint induced by the short-range of the potential (a > 18) facilitates a progressive 
layering. For example, for the case iV = 13, the icosahedron structure, which is particularly 
stable in the LJ case [lflj . is not observed in the short range case. 

It is instructive to compare these results also with hard spheres clusters geometries, 
calculated theoretically minimizing the second moment of the mass distribution M [20, 25] , 
and recently measured with a new experimental technique introduced to produce compact 



clusters of controlled numbers of small colloidal particles J26 
have been found to be identical to those calculated in Ref. 



For N < 11 these clusters 
20j. The third row of Fig. El 



reports the HS clusters from Ref. 20]. As for the comparison with the LJ case, for N < 7 
all clusters have the same structure, while for N > 7 differences in the potential start to be 
significant. 



III. THE COMPETITION BETWEEN ATTRACTION AND REPULSION 
TERMS 

In this section we report calculations of the cluster ground state energy for a potential 
composed by a short-range attractive part (modeled only with a = 100, since we have shown 
previously that a = 18 can not be distinguished from the a = 100 case) complemented by 
a screened electrostatic repulsive interaction modeled by a Yukawa potential, 

Mr) = A ^- (2) 

The resulting total potential is thus V(r) = V2a-a +Vy. In the following, A is given in 
unit of e and £ in unit of a. Fig. |3] shows the shape of the total potential for different 
£ and A values. The resulting potential is conceptually similar to the well-known DLVO 
potential J23| , even if it differs for the absence of the weak secondary minima and the finite 
value of the attractive interaction energy. 

Fig. 13] shows the ground state energy per particle E/N for different A and £ values. The 
same minimization procedure as described in the previous section has been implemented, 
with the additional condition that, in the search procedure, moves creating disconnected 
clusters are rejected. For appropriate values of A and £ (for example A = 0.2 and £ = 2.0), 
the cluster ground state energy shows a minimum at a finite size N*, indicating that clusters 
of size larger than N* are energetically disfavored. A simple physical explanation of the 
existence of an optimal size can be given as follows. At large enough N, the addition of 
an extra layer of particles will contribute to the energy with a negative term due to the 
attractive nearest neighbor interactions (which involves only interactions with particles in 
the surface layer due to the short range of the potential), plus a positive contribution arising 
from the Yukawa repulsion which, due to its longer range, involves instead a large fraction 
of particles in the cluster. The balance between these two terms provides a condition for 
iV*. Of course, the larger the amplitude or the longer the range of the repulsion, the smaller 
iV* is. The fact that an optimal finite cluster size is observed, suggests that a macroscopic 
system at low temperature will not form a single aggregate, but will prefer to partition the 
particles into clusters of size N*. In this respect, liquid condensation is inhibited and the 
structure of the system at low temperature and low packing fractions will be constituted by 
a cluster phase. At finite temperatures, entropic contributions to the free energy will become 



important and will always favor the stability of clusters of size smaller than N*. Hence, N* 
plays the role of an upper limit for the cluster size. In other cases, a minimum is not found, 
but the energy per particle becomes essentially flat. In these conditions, the ground state 
of a macroscopic system will be composed by a highly polydisperse cluster phase. 

In Fig. we compare the ground state cluster geometries for N = 38 calculated for dif- 
ferent choices of the Yukawa parameters. At fixed cluster size, by changing the potential pa- 
rameters the cluster structure can be tuned from almost spherical to almost 1— dimensional. 
The studied parameters encompass situations occurring in uncharged colloids (A = 0) with 
those typical of charged colloids (A > 0) in apolar solvents or weakly screened polar sol- 
vents. Indeed, A is proportional to the effective charge of the colloidal particle, while the 
screening length £ is controlled by the ionic strength. In the absence of salt, £ depends 
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loid concentration, and a small screening is produced by the counterions in the 
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29]. 

Fig. El shows the iV-dependence of the cluster geometry, for the case A = 0.05 and £ = 2, 
to highlight the changes in the preferential geometry with cluster size. It appears that the 
geometry of the clusters is strongly size dependent. To convey this point, we show in Fig. 
the size dependence of the gyration radius Rg, defined in terms of particle's coordinates r, 
and center of mass coordinates Hcm as 
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For spherical clusters, Rq ~ iV 1 / 3 , for planar structures, Rg ~ iV 1 / 2 , and for linear struc- 
tures, Rq ~ N. Fig. [7| shows that, while in the pure short-range attractive case clusters 
retain their spherical shape for any size, in the case of addition of a repulsive potential only 
clusters of small enough size are spherical. On increasing N, the ground state cluster struc- 
ture becomes more and more linear. The sharp cross-over size between spherical and linear 
cluster shape decreases on increasing A or £. Similarly, the linear cluster becomes thinner 
and thinner when repulsive effects becomes more and more relevant, as evidentiated by the 
amplitude of the linear iV-dependence of Rq- Once the linear cluster shape is established, 
the cluster energy increases linearly with the the cluster size, resulting in the essentially flat 
iV- dependence of the energy per particle (see Fig. HJ). 

Three possible scenarios appear to take place in the cluster structure of the studied 
potential: i) very small A values, where the repulsive energy is not sufficient to overcome 



the attractive contribution. Under these conditions, clusters are spherical and the standard 
behavior (infinite optimal size, liquid phase) is recovered; ii) an intermediate case, where a 
minimum in the cluster energy per particle is observed at an optimal size N*, with a spherical 
cluster structure (as, for example, the case A = 0.2, £ = 2). Linear clusters can be built 
but with an energy slightly higher than that of the optimal cluster size; in) a monotonically 
decreasing cluster energy per particle, ending into a flat curve, signaling equivalent stability 
for linear clusters of very different size (as, for example, the extreme case A = 7.93, £ = 0.5). 
It is tempting to speculate that when parameters are chosen in such a way that linear 
growth is preferential, at finite but small temperature (i.e., when structures different from 
the ground state structure are also probed) clusters made of arms branching out from regions 
of locally higher energy can be generated. These structures, if macroscopic in size, would 
generate a gel-like structure, since no driving force for macroscopic aggregation is present 
(due to its inhibition caused by the repulsive term). 

IV. THE SPHERICAL CASE: THE CROSS-OVER FROM INFINITE TO FINITE 
OPTIMAL SIZE 

The numerical results presented in the previous sections confirm the possibility of cre- 
ating stable clusters of optimal size, playing with the competition between the short range 
attraction and the long range repulsion. We have discussed the fact that in the limit of low 
temperatures and small packing fractions, i.e., in the conditions where entropic contribu- 
tions as well as cluster-cluster interactions can be neglected, the cluster ground state energy 
E plays the role of relevant thermodynamic potential and, therefore, the minimum in E/N 
vs N provides an estimate of the stable cluster size N*. Under these conditions, the system 
will partition into clusters of size N* and no thermodynamic (macroscopic) coagulation will 
take place. Hence, the condition N* = oo acts as critical condition for the existence of a 
macroscopic phase separation as opposed to a phase made of clusters of finite size N*. 

To estimate the critical line separating the liquid phase (iV* = oo) from a cluster phase 
(iV* finite) in the (A, £) parameter space, it is necessary to calculate the ground state energy 
for clusters of large size, which prevents the use of the "exact" numerical minimization which 
have been presented previously. An estimate of the critical line in the (A, £) space can 
be calculated analytically, under the hypothesis that clusters have a spherical shape and a 



homogeneous density of particles. The assumption of a spherical shape is expected to be valid 
in the region of parameter space where N* = oo, i.e., where A or £ are small and attraction is 
still dominant, as shown in the previous section. The numerical study of the iV- dependence 
of the ground state of the 2a — a potential reported above shows that the iV-dependence 
of the attractive part of the potential can be written as V2 a -a{N)/N = Cq + C\N~ 1 ^ with 
c ~ —6 and c\ ~ 7.4. Since the value of the bulk energy (co) reflects the energy of a very 
compact state (in which each particles is surrounded by 12 neighbors), the local density can 
be approximated with the local density of the fee structure, providing a mean to convert the 
cluster radius R to the cluster size using N = r]f cc (2R/a) 3 , with r]f cc ~ 0.74. 

To estimate the repulsive energy of a cluster of radius R 3> er, we assume for simplicity 
a particle-particle radial distribution function g(r) = 0(r — cr) where is the Heaviside 
step function. This choice ensures that unphysical interactions, which would arise by pairs 
of points in the sphere closer than a, are eliminated. A more precise calculation could 
be performed numerically by implementing a more detailed expression for g(r) (which, for 
example, could be calculated with the Verlet- Weiss g(r) for HS, or with the known g(r) of 
the fee structure). The simpler step- wise approximation used here is sufficient to predict 
the shape of the critical line. For large cluster sizes, the resulting expression for the cluster 
repulsive energy per particle is (see Appendix ^ 

-24e(R + 2 , 16(cr + i2(a 2 + 2< + 2£ 2 ) 
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Fig. |H1 shows the 1/R dependence of the cluster energy per particle, i.e., the sum of the 
attractive and the repulsive part, for £ = 1 and several values of A. On increasing A, a 
progressive bending up of the curves takes place for large cluster radii, until for a critical 
value A c (ss 0.17 for £ = 1), the slope of E(R) becomes flat at R — oo, signaling that 
the lowest ground state does not require that all particles belong to the same cluster. For 
A > A c , a minimum of E(R) arises at finite R value. To estimate the critical value of A(£), 
the expression Eq. (jlj) for Ey/N can be expanded in powers of 1/R. At first order in 1/R 
we find, 

Ey(R) A A £ 3 (a \ a A £ 4 fa 2 a \ * 1 



N 



^ A Vfcc^[^ + iy^-l8A Vfcc ^^- + 2- + 2Je-T-, (5) 



which provides, by adding the attractive part (also linear in 1/R) and setting the resulting 
coefficient of the 1/R part to zero, A c as a function of £, 

c a 4 £ 2 

Ac ^ = 36^ e W + 2£a + a 2 ^ ' (6) 

The critical line A c (£)) shown in Fig. El in the (A,£) plane, locates the region where, at low 
temperature, a phase of clusters of finite size is expected. For £ > a, A c ss £~ 4 . 

Finally, to provide an estimate of the region where, even at T = monomers are the stable 
state, we show in Fig. 01 the line corresponding to the condition of the energy of a dimer Ed 
being zero, i.e., along this line the pair potential repulsive energy at distance a is equal to the 
short-range attractive energy (i.e., Ae~ a ^/(cr/£) = 1). Crossing the E^ = line from below, 
the dimer ground state energy goes from negative to positive values. The intersection of this 
line with the A c (^) line suggests that, for £ < 0.2, there is no possibility of a cluster phase, 
with clusters larger than monomers. Therefore, for such small values of £, the ground state is 
given by an infinite cluster when the dimer energy is negative, crossing sharply to a solution 
of monomers when the dimer energy becomes slightly repulsive. Physical realization of 
cluster phases should then be searched in systems where the screening length is comparable 
to the particle size, i.e., in weakly polar solvents or for small colloidal particles. 

In concluding this section, it is worth recalling that, as seen in the previous section, 
for sizes comparable or larger than N*, the shape of the cluster becomes essentially linear 
and the spherical cluster calculations lose progressively their meaning. Hence, the spherical 
calculations reported here should be limited to the case iV* = oo, i.e., below the critical 
Ai(0 line. Close to the critical line and above, more refined calculations accounting for 
almost linear cluster shapes should be performed. 

V. CONCLUSIONS 

In this paper, we have studied, for a model where short-range attraction and long-range 
repulsion are simultaneously present, the structure and the energy of the ground state for 
clusters of different sizes. In these systems, microphase separated states --in the form of 
spherical (micellar) or columnar phases — arise spontaneously, due to the introduction of a 
new characteristic length provided by the balance between attractive and repulsive energies. 
Under the hypothesis of spherical clusters, an analytic evaluation of the cluster energy has 
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been performed, providing a criterium for the existence of a cluster phase. 

By varying the two parameters of the repulsive Yukawa potential, controlling respectively 
the amplitude of the repulsion and the screening length, clusters of different morphology, 
from almost spherical — for small A and £ values — to almost one-dimensional — for larger 
A and £ values — can be generated. The evaluation of the ground state energy as a function 
of the cluster size gives evidence of a progressive tendency toward one-dimensional growth 
for all cases leading to cluster phases. This preferential one- dimensional growth is expected 
to enhance the stability of collective ordering into columnar or lamellar phases when cluster- 
cluster interactions are taken into account, both energetically and entropicallv 30L in full 
a„t With the predict for unsC reeoed repulsive inactions fl MM}. 



It is tempting to connect the one- dimensional growth followed by a dynamical arrest 
phenomena, which is observed in several protein solution systems J3l|, |3j, |33|, |34|, to the 



results discussed in the present article. Indeed, for these protein solutions, a change in 
temperature or in the solvent properties can trigger an aggregation process of proteins 
into cylindric clusters. Under appropriate concentration conditions, these one-dimensional 
clusters further associate to form a macroscopic gel. The mechanisms discussed in this 
manuscript account for both the formation of cylindrical clusters and the insurgence of an 
effective vanishing surface tension, a condition necessary to stabilize a gel-phase with respect 
to phase separation. 

Finally, we want to stress that in this study we have focused on isolated cluster properties, 
and attempted to connect the cluster properties to the formation of a cluster phase, as 
opposed to a condensation of a dense liquid. In doing so, we have neglected the cluster- 
cluster interactions, which will play a very relevant role also at low temperatures, for packing 
fractions at which the cluster-cluster distance becomes comparable with the screening length. 
This will bring into play not only thermodynamic considerations but also, due to the low 
temperature, kinetic considerations. Preliminary work in this direction indeed suggests that 
slow-dynamics phenomena, related to the cluster-cluster repulsive interactions, may play a 
relevant role in arresting the equilibration of systems of particles interacting with the type 
of potential studied in this work [15]. 



11 



Acknowledgments 

Authors acknowledge support from MIUR Cofin 2002, FIRB and MRTN-CT-2003- 
504712. We thank A. Scala and J. Groenewold for helpful discussions. 



12 



APPENDIX: CALCULATION OF THE YUKAWA CLUSTER ENERGY 



To evaluate the cluster energy, under the hypothesis of a homogeneous spherical cluster 
of number density p, we proceed in two steps: First we evaluate the potential energy of a 
particle located at distance x from the center of the sphere of radius R, by integrating all 
contributions from points located at the same distance r from the selected particle. Second, 
we evaluate the cluster energy by summing over all particles located at distance x, integrating 
x from to R. The energy of a particle located in x inside a sphere of radius R, W(R,x), 
can be written as 



W(R,x) =p 



R—x pR+x 

An I drr 2 Vy(r)g(r) + / drS(r; R, x)Vy(r)g(r) 

JR-x 



(A.l) 



where S(r; R, x) is the surface generated by the intersection of two spheres of radius r and 
R, whose centers are located at distance x apart. From standard geometry, 

- 2 - r 2 + R 2 



S(r; R, x) = 2irr 



x 



2x 



(x — r) 



(A.2) 



The condition g(r) = 0(r — a) acts in the integration limits providing two solutions, one for 
the case R — x < o and one for the case R — x > a, 



-R+x 



p J a x drS(r; R,x)Vy{t) 



x > R — o 



W{R,x) = < 



(A.3) 



p An f*~ x drr 2 V Y {r) + f^ drS(r; R, x)V Y (r) 
The total potential energy of the cluster Ey{R) is thus calculated as 

f R 

E Y {R) = Anp / dxx 2 W{R,x). 
Jo 



x < R — a 



(A.4) 



For the case of the Yukawa potential, standard integration provides the result for Ey(R) 
reported in Eq. (J1J). Indeed, after converting to energy per particle and using |cr 3 p = r]f cc , 
one finds, 



E Y {R) 3A Vfc 2 



N 



2a 3 
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FIG. 1: Ground state energy per particle as a function of 1/iV 1 ' 3 (N on the top x-axis) for the 
2a — a potential at the indicated values of a. The values of the energy per particle for the fee 
crystal structure are indicated by arrows. Inset: Radial dependence of interaction potential V2 a -a 
for the three studied values of a. 
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FIG. 4: Ground state energies per particle for five different choices of the potential parameters. 
The energies have been calculated by basin- hopping Monte Carlo optimization []j|, as discussed 
in the text. 
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FIG. 5: Ground state clusters for N = 38 for the indicated potentials. On changing the values of 
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state structure. 
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FIG. 6: Ground state clusters for TV = 4, 8, 16, 32, 64 for A = 0.05 and £ = 2.0. 
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